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In recent years we have witnessed a real increase in the production of wind 
turbines and wind farm installations around the world, in order to improve 
this own energy, several studies have focused on the interest of controlling 
the active and reactive power of the system. Wind power, and at the same 
time on the quality of the energy produced and its connection in order to 
ingest suitable electrical energy into the distribution network. This article 
studies a new control technology to meet the various constraints in the field. 
The objective for work is to develop and study the sliding mode control 
method applied to a wind power system based on doubly fed induction 
generation (DFIG), as well as an optimization using the fuzzy logic 
technique. Ensuring the stability of the system is one of the objectives of 
using the Lyapunov nonlinear technique in the sliding mode control strategy 
which will be applied to the two converters (machine side and network side). 
The proposed solution was to validate a simulation on MATLAB/Simulink, 
tracking test (true wind speed) and also the robustness of the system. 
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NOMENCLATURE 
WPGS Wind power generation system Tar, Tgr Rotor Current 
DFIG : Doubly fed induction generator as, Pgs Stator flux 
Pr : Power ricovred by the wind turbine ar, Oar Rotor Flux 
Pyent : Kinetic wind power Ps : Active Stator Power 
Cp : Power coefficient Qs : Reactive Stator Power 
Tr : Aerodynamic torque Pr : Active Rotor Power 
Vas, Vas Stator voltage Qr : Reactive Rotor Power 
Var, Var Rotor voltage Tem : Electromechanical torque 
Tas, Tos : Stator current 


1. INTRODUCTION 


In recent decades, the massive industrialization caused the hard increase of electrical energy 
consumption. The forecasts of energy needs for the coming years only confirm this. Finding different 
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sustainable solutions is necessary to meet the future word energy needs [1]-[3]. In the immediate future, the 
world experiences a recourse to renewable energies, which present inexhaustible resources. Among which we 
find wind energy identified as one of the most promising for achieving greenhouse gas emission reduction 
targets. The exploitable renewable energies are numerous and varied. That said, the development of wind 
energy is booming, so wind systems are developing more quickly because we can now rely on reliable and 
proven techniques with better performance. 

Wind energy therefore remains the most used in energy production technique, hence the case of 
Morocco, which relies on wind energy production as the primary source of renewable energy with 4699.4 GWh 
in 2019 against 1581, 50 GWh for solar and it is considered the country's second energy source against 
32,218.50 for thermal energy, i.e. a rate of 11.7% of the country's total energy production [3], [4]. The Figure 1 
shows the graphic of energy production in Morocco. 

Currently, the majority of wind turbine installations in the world used DFIG- Generator. Because 
this type of machine has the ability to produce electric power at variable speed. The control of the studied 
system based on double-fed machines by classical methods (proportional, integral and derivative action) 
implies the knowledge of the various parameters, which can influence it and allow the system to achieve the 
set results. In addition, the strongly coupled variables, increasing the system, makes control more difficult. 
Using nonlinear methods is one of the solution to solve this problem like the sliding mode controller (SMC) 
control. The latter is considered one of the simplest approaches to control nonlinear systems and systems 
with imprecise models [5]—[10]. 
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Figure 1. Graphical representation of energy production in Morocco 


2. WIND ENERGY CONVERSION SYSTEMS (WECS) MODEL 
The Figure 2 presents the wind energy structure. It contains the doubly fed induction generator, the 
converter, the filter and the grid. The control is showing in the input of the converter. 
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Figure 2. WECS structure 


2.1. Resource modeling 
2.1.1. Resources potential 
The equation for the wind power recoverable by the wind turbine is [10]—[13]: 


Pr = Prent: Cp (A,B) (1) 
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The total kinetic power in the wind turbine is: 
22 3 
Prent = 2 pSv (2) 
Pk = 5S. Cp (A, B)v3 (3) 


For the production of wind turbines, it is necessary to take into account the defined range of the energy 
extraction coefficient which is generally between 0.35 and 0.5. 


1 =C 1 
CA B) = cy. (c24 - c3. B- c4) -€ A+ C,¢.A (4) 
1 1 0.035 
A At+0.08B 1+6? (5) 
— WR 
A=" (6) 
The equation of the aerodynamic torque (Tr) whose coefficient curves. 
Py 1 1 

ee ae z PS. CoA, Pv- (7) 


The following two equations used to model the adaptation of the two turbine speeds to that of the generator. 


0=G.0, 

T, = G. Taq 

doa. 

JË = Tm- Tem — fN (8) 


It is necessary to have a simulation tool capable of modeling the entire energy conversion chain and 
predicting its performance. Initially, a model of the wind generator is proposed under MATLAB/Simulink to 
evaluate its dynamics and its performance under different operating parameter changes. Figure 3 shows the 
model of bloc turbine in MATLAB/Simulink. 
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Figure 3. Turbine model 


2.2. Pitch control and mechanical control 
In order to better control the operation of the machine, a number of parameters must be defined 
beforehand on which the operation and behavior of the machine depends, such as: 
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— The wind speed which allows defining the parametric changes of the machine, 

— The nominal power P,on the wind turbine shaft, 

— The speed reference v, which corresponds to the power P, 

— The starting speed vg is the minimum wind speed can start to operate and provide energy. 

— The safety speed Vm is the wind speed which when it is exceeded the wind turbine automatically stops 
working. 

The turbine operating power diagramed is divided into four zones defined according to the wind 
speed Vn, Vd and Vm witch presented in Figure 4 the first zone corresponds to P=0 the turbine is not yet 
working. The second zone corresponds to the phase where the power supplied depends on the wind speed. 
The third is when the supplied power stabilizes at Pn for a constant speed of rotation and in the end the fourth 
zone which corresponds to the stopping of the energy transfer for safety reasons [14]-[19] 

Figure 5 shown the turbine has greater efficiency for low stall angle values. A single control loop may 
be sufficient to control the pitch for angle B to manage the orientation of the wind turbine blades. The rated 
power is controlled by a single servo loop and which is read with an open loop for controlling the orientation of 
the blades or the output of the first is the values of the pitch angle of the second. and therefore, we have two 
servo loops, one for nominal power control the second is for wedging angle control. 


0 Va Va Vm Vms) 


Figure 4. The turbine operating power diagramed Figure 5. Pitch angle influence 


2.3. Power control technic 

In order to limit the extracted power from the turbine so that it is equal to the nominal power of the 
generator, two principles of aerodynamic control can be used [20]-[24]: he first is the "pitch control" or not 
variable "which is based on the rotation of the blades by electric or hydraulic actuators which makes it possible 
to adjust the pitch angle B, which is called the blade orientation mechanism. it is therefore a principle which 
allows the mechanical adjustment of the turbine in order to keep it at its maximum efficiency. 

The second principle is the "aerodynamic stall" considered more robust and the simples. The 
moment when the wind speed exceeds the nominal speed of rotation of the turbine, a suction thus created by 
the inclination of the blades which will not be in the same optimal position leads to a reduction in the power 
coefficient. Therefore, this inclination limits the increase in the speed of rotation. The interest of the first 
check can be observed from Figure 6. 
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Figure 6. Power of a turbine 
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2.4. Generator modeling 

The generator we used for wind power conversion is the dual-feed asynchronous generator DFIG. 
[25]-[27]. The DFIG Model is schematized in Figure 7 is a in the park referential model, using the following 
equations: 


Couple électromagn 


Pr 


eqs magnétique 


courants de sortie 


Figure 7. DFIG model 


2.4.1. Stator and rotor voltage 
The following equations define respectively. The stator and rotor voltage related to stator and rotor 
flux. We use in this definition park referential. 


das 

Vas = Rslas + Eds = Os Pqs 
d@qs 

Vas = Rslqs + — 0s Pas (9) 
d r 

Var = Rrlar + far = 0-Pqr 
d@ qr 

Var = RyIgr + a — 0,@ar 

2.4.2. Stator and rotor currents 


Using the park referential. We define the following equations respectively. The stator and rotor 
current related to stator and rotor flux 


I 1 M 
ds 7 OL Pas oLy Par 
1 M 
las = i = 10 
as T SL; Pqs ols, Pqr (10) 


lar = —. Qar- —. 
dr o.Ly Gar o.Ly Pas 


I . —_-—. 
ar 5 g Pqr SLL Pqs 


2.4.3. Stators and rotor flux 
In the same referential dq we will define. The stator flus related to stator and rotor currents. The 
rotor flux related to stator and rotor currents and the electromagnetic torque relative to rotation speed 


Qas = Lslas + Mlar 
Pesk ket Mi (11) 
Qar = Lylar + Mlas 
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gr = Lylgr + MIgs 

Tem = Ty + £0 + JS (12) 
2.4.4. Stator and rotor powers 

The (13) describe the active and reactive power of the stator in a two-phase reference system 


relative to current and voltage. The active and reactive power of the rotor in a two-phase reference system 
relative to current and voltage. 


P, = = (Vas- las + Vas- las 
Qs = $ (Vas: las + Vas- las (13) 
P, = Ž (Var- lar + Vor lqr 
Qe = 5 (Var-lar + Var- lar 


2.5.2. DFIG model 
The following equations determine the electromechanical torque Figure 7. 


Tem = Pas- Pqs — Igs- Pas) (14) 


The reference frame used is the d-q reference frame related to the rotating field of the stator, and by 
considering the hypothesis for the generators used in wind energy the resistance of the stator will be 
neglected, this said the equation of the torque will be: 


M 
Tem = PI (Pas: Igr = Pqslar) (15) 


3. SLIDING-MODE-CONTROL APPLIED TO THE DFIG 
This method consists in first determining the adequate sliding surface, as well as the equivalent and 
non-linear value for each quantity to be regulated [27]-[29]. 


3.1. Sliding surfaces 
x = f(x, t) + B(x, t). u(x, t) (16) 


Is a representation of a non lineair système 
xeR” 
f(x, t)eR", B(x, t)eR™™, ueR™ 


f and g are two continuous nonlinear functions uncertain and bounded Considering the "J. J Slotine" 
equation, the sliding surface [28] is: 


s(x) = (+8) ec) 17) 


where: 

e=Xq-x 

6 : Positive coefficient; 

e : Error to be regulated; 

n: System order. 

Xq : Desired variable 

This command is based on the principle to keep the surface at zero, so we can have just one solution of the 
differential equation e(x) = 0 . 


3.2. Conditions of convergence 
The convergence condition is defined by the lyapunov equation [29]. 


V(x) = 5.S()? > 0 (18) 
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3.3. Determination of the control law 


The command (u) is a variable structure command and the control algorithm is defined by the 
relation [29]. 


= es s(x,t) > 0 
~ tu~ (x), s(x, t) <0 
U = Ueg + Un (19) 


Where: 
u` + ut Where utand wu: continuous functions 
The derivative of the surface S (x) is: 


S(x) = s, f(x) + s B(x). Ueq + A B(x). Un (20) 
-1 
Ueq = -% 4@).(%.B09) (21) 


as 
So =. B(x) #0, 
By in S(x): 
: as 
S(x) = z BO): Un 


Considering the attractiveness condition S(x).S(x) < 0 , the sign U,, will be opposed to H B(x). S(x) 


Different form existed in the literature [28], like U, = K.sign(S(x)) (22) 


K : Positive constant 
(sign): the classical sign function. 
The saturation function can be expressed by [27], [28]: 


lifo>e 
sat(”) y ifp <e (23) 
zif lol se 


3.4. The sliding mode application 
This part is devoted to the application of the control by sliding mode in the asynchronous generator 
with double feed. Following the new considerations, the previous become as such: 


Vsq = 0 


Vsq = Vs = Ws. Øs 
do 
Vrda = R,. La + ae + Ory. Ø 


rq 

Veq = Rr- Irq + et + 00r- Ora (24) 
Bei ha (25) 
Qs = $ Veira (26) 
İra = = = a + Wr Irq (27) 
İrq = = p a [rq 7 Or. moo (28) 
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4. ACTIVE POWER CONTROL SURFACE 
Considering the sliding surface by J. Slotine given by: 


s(x) = (4+8) ec) (29) 


We consider n = 1, 
We have the active power control surface the: 


S(P;) = e(P;) = Pyrep — Ps 
S(Ps) = Poret — É, (30) 
We replace P, by its expression as well as ley, we get: 


MV, 


. R M A 
S(P,) = Pref + Vs- Le ee Lro’ Irq Wr. LyLpJ.05 (31) 
7 * M Vr e +Vr n Rr M. Vs 
S(P;) = Pref + Vs- i, sa = a Le Teg — Or- Mere (32) 
In sliding mode: 
S(P;) = 0 
S(Ps) =0 
Vign = 0 
. Vs 
Vige Le Lea" Poret + (Rr Irq + Lr: 0. Ira + @r- M. =) (33) 
The control is given V-qn by: 
Vrqn = Ks. sat(S(Ps) (34) 
With: K,: Positive gain 
5. REACTIVE POWER CONTROL SURFACE 
With the same method as the active power, we find: 
S(Qs) = e(Qs) = Qsref =Q; 
S(Qs) = Qsref — Qs (35) 
We replace Q,by its expression as well as Iq, we get: 
: . M Vr Rr 
S(Qs) = Qsref + Vs- Le’ = SRo Tra (36) 
By replacing the expression of V,.qwith V,geq + V,qy, we obtain: 
z ; M Vrdeq+tVrdn Rr 
S(Qs) = Qsref + Vs. Le seas eas = Laird = Ory. Irq (37) 
In sliding mode: 
S(Qs) =0 
S(Q;) =0 
Vian = 0 
V,deq becomes: 
Videu = —Ls. Le- av Qsref + Rr: Ira — Lr- O. Teg: Or (38) 
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and discontinuous control V„anis given by: 


Vran = Kg. sat(S(Qs) (39) 


with: Ka: Positive gain 


6. SLIDING FUZZY MOD CONTROL 

The sliding mode control maintains robustness and improves system performance. However, the 
command switching term produces a chattering phenomenon which is its main flaw. overcome this problem, 
we propose the model that combines between fuzzy logic control "FLC" and sliding mode control "SMC" 
[30]. The sliding mode control law is defined by the control strategy, which depends on the sign of the error 
if it is negative, and then the system output is driven in the positive direction. 

For this the term Kg.sat(S(U,) can be replaced by a fuzzy controller. This controller has an input 
S(x) and an output U, and its rule base is used to establish a connection between S(x) and Up. This is 
interpreted by the rules of the form (If... Then) [31], [32]. We design for this purpose a command by fuzzy 
sliding mode to control the active and reactive powers as presented in Figure 8 The basic input and output 
linguistic variables of the fuzzy control rule are presented in Figure 9. 


Figure 8. Membership functions defined for sliding fuzzy control of reactive power 


Commande par mode glissant flou 


Commande 
floue 


Figure 9. Active power control principal diagram by fuzzy sliding mode 
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In this section, the FLC fuzzy command is introduced to replace the Kg.sat(S(U,) function, good 
steady-state dynamics can be obtained applying fuzzy logic to the DFIG speed command by combining it with 
the sliding mode command. The verification of the advantages of this new fuzzy sliding mode command will be 
carried out by simulation. The speed and rotor flux regulators are replaced by a sliding-fuzzy mode regulator. 


7. SIMULATION RESULTS 

In purpose to present the evaluation of the performance of the new hybrid fuzzy sliding mode 
control SFMC, we will expose the simulations carried out on an asynchronous machine with double power 
supply, Figure 10 represents the block diagram of the simulation of the principle of the fuzzy sliding mode 
control applied to DFIG, replacing the sign (Sw) function with a crazy controller. 


Parameter 


P, (2) Q 6) estimation 
P, set Qu cet 


Active and reactive 
power controllers 


Reference power on 
rotor side 


Figure 10. Sliding mode control applied to DFIG 


7.1. Performance test 

We will study and analyze the response of the DFIG for a wind profile that we choose to be random. 
In order to improve the quality of the electrical energy injected into the network, we are going to set the 
reactized power to zero. (Qs_ref = 0 Var). The wind profile used in the simulations presented is a real wind 
profile. The results of the responses of the fuzzy logic mode sliding control of the DFIG are represented in 
Figure 11 (see Appendix) shows the chosen winds are that of the Tanger region located in the north of 
Morocco. The Figure 12 shows THD for the rotor current 6.29% and 0.16% for the stator current 
(according to standard IEEES19-2014). 


i a (50Hz) = 9.903 , THD= 0.16% 
T T 


ig (Hz) = 4.05 , THD= 6.29% 
r H 5p r : H 


Mag (% Fondamental) 


H H i ; ò H : 
0 400 800 1200 1600 2000 0 400 800 1200 1600 2000 
Frequency (Hz) Fréquence (Hz) 


Figure 12. THD of currents at variable wind speed 
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7.2. Robustness test 

We are going to make some modifications to the internal parameters of the DFIG to verify the 
robustness and stability of the fuzzy sliding mode control on active and reactive power regulation. On this, 
the tests carried out are carried out on the basis of an increase and decrease of the parameters of the machine 
such as the stator, the resistances and the inductances of the rotor. Thus, the dynamic behavior of the system 
is represented in Figure 13. 
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Figure 13. Robustness test during resistance variation 


However, the changes applied to the resistances of the stator and the rotor cause a slight increase in 
the response time, which is detected in Figure 13, we therefore notice that the simulation results show a great 
robustness of the fuzzy-sliding mode command with respect to the parametric variations of the DFIG. This is 
justified by the good follow-up of the power setpoints, with almost the same response time at start-up. It 
always maintains the decoupling between active and reactive power. This performance is showing in Table 1. 


Table 1. Performance comparison 


Publication Technic Control Efficiency Error Coso « Robustness » 
[10] Adaptive Controller 93.5 0,15 % 0.997 Moderate-high 
Fuzzy Controller 93.99 0,14 % 0.974 Moderate-high 
[4] classical DTC 92.13 0,32 % 0.983 Moderate-high 
PI-DTC 92.07 0,12 % 0.978 Moderate-high 
[3] Sliding mode Control 94.82 0,2 % 0.972 Low 
Proposal technique Fuzzy-SMC Control 98.99 0,12 % 0.995 High 
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8. CONCLUSION 

In this article, we have presented the control of a wind energy recovery system equipped with a 
doubly fed induction generator. Firstly, a generator model was proposed, and then a fuzzy sliding mode 
control strategy of the asynchronous generator allowing independent power control was proposed. Fuzzy 
sliding mode active and reactive power regulators have been proposed and tested. The simulation results 
enabled us to assess the performance and robustness of the fuzzy sliding mode control. Thanks to the 
response characteristics, good performance is observed even in the presence of set point variations. The 
pursuit of power is without overshoot. Decoupling, stability and convergence towards equilibrium are 
ensured. In addition, this setting has the advantage of being easily implemented in a computer control. 


APPENDIX 
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Figure 11. Performance test of the SFMC applied to the DFIG 
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